arXiv:1503.02407v2 [astro-ph.CO] 7 Apr 2015 


ScalPy: A Python Package For Late Time Scalar Field Cosmology 


Sumit KumaiQ and Anjan A. Sei0 
Centre for Theoretical Physics 
Jamia Millia Islamia, New Delhi 110025, India 

Abhishek Jang@ 

Department of Physical Sciences 
Indian Institute of Science Education And Research Kolkata 
Mohanpur, West Bengal 7fl2f6, India 

We present a python package ’’ScalPy” for studying the late time scalar field cosmology for a 
wide variety of scalar field models, namely the quintessence, tachyon and Galileon model. The 
package solves the autonomous system of equations for power law and exponential potential. But 
it can be easily generalized to add more complicated potential. For completeness, we also include 
the standard parameterization for dark energy models, e.g. the ACDM, wCDM, rooWaCDM as 
well as the GCG parameterization. The package also solves the linear growth equation for matter 
perturbations on sub-horizon scales. All the important observables related to background universe 
as well as to the perturbed universe, e.g. luminosity distance ( Dl{z )), angular diameter distance 
( Da{z )), normalized Hubble parameter ( h(z )), lookback time (U), equation of state for the dark 
energy ( w(z )), growth rate (/ = ^ s a ), linear matter power spectra (P(fe)), and its normalization as 
can be obtained from this package. The code is further integrated with the publicly available MCMC 
hammer “emcee” to constrain the different models using the presently available observational data. 
The code is available online at https://github.com/sum33it/scalpy 


I. INTRODUCTION 

There is compelling evidence that the universe is currently undergoing a phase of accelerated expansion [lH3| • This 
indicates some new physics at cosmological scales (“dark energy”) [§J that can give rise to repulsive gravity resulting the 
universe to accelerate. The cosmological constant with w = ^ = — 1 is the simplest example of repulsive gravity and is 
consistent with the majority of cosmological observations. But the observational data at present is also consistent with 
dark energy behavior which is strictly not constant but evolves slowly throughout the cosmological history. Moreover 
the cosmological constant comes with its own baggage of problems like cosmic coincidence problem which can be 
avoided using dynamical dark energy. 

Scalar fields are the most widely used dynamical dark energy models where late time acceleration can be obtained by 
adjusting the slope of the scalar field potential around suitable epoch [S], But the cosmological evolution of these models 
is severely constrained by very accurate cosmological observations. From the measurements of temperature anisotropy 
in the cosmic microwave background (CMB) radiation, the distance to last scattering is very well determined. This 
restricts the equation of state for the scalar field to be very close to w = — 1 at present. The type-la supernovae 
observations as well as Baryon Acoustic Oscillation (BAO) measurement in the large scale structures also give similar 
constraint on the scalar field equation of state. On the other hand CMB measurements also restricts the scalar field 
contribution to the total energy density of the universe during recombination to be less than 1%. All these constraints 
disfavor 0 the original “tracker” kind of quintessence model [ll] where the scalar field initially has a fast roll phase 
before settling to the cosmological constant behavior at present. Just recently, analysis of the spectrum from a distant 
quasar finds no evidence of deviation in molecular lines produced around 12 billion years ago. This confirms no 
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change in the mass ratio of the proton to the electron. The results also indicate that a dark energy scalar field— if it 
exists- has not evolved appreciably over 90% the age of the Universe jl2j], thereby putting doubt on the tracker kind 
of models as a possible candidate for dark energy. 

The other scalar field behavior that is consistent with the present observational data is the “thawing” one [Tl| . 
Currently there are many types of scalar field models that can give rise to thawing behavior. Examples are quintessence 
tachyons [l3j as well as Galileon [14] models. In all these thawing models, the thawing nature does not depend 
crucially on the form of the potential. Once the initial conditions are chosen in such way that the scalar field starts 
in a nearly frozen state with w ~ —1, the thawing behavior is ensured. Fine tuning this initial condition is an issue 
in all thawing type of models. In a recent investigation, attempt has been made to relate this initial condition to 
inflationary dynamics |l(i. 

Given the fact that thawing models are observationally promising as a possible candidate for dynamical dark energy 
models, we need to have a combined setup to investigate the cosmological evolution for different thawing scalar field 
models. In this work, we present such a setup in a Python environment. We name this setup as “ScalPy” which 
solves the autonomous system with a given scalar field model by providing suitable initial conditions and calculates 
the observables like normalized Hubble parameter, angular diameter distance, luminosity distance, growth function, 
growth rate, linear matter power spectrum etc. We consider two simple potentials , the power-law and exponential 
for all the different versions of the scalar fields as these are the two most widely used potentials for dark energy in the 
literature. But one can easily incorporate the other potentials in this package. We also combine the package with the 
publicly available MCMC sampler “emcee” E2 to do the statistical analysis of different scalar fields with observational 
data. 

Given the promise of thawing scalar field models as possible candidate for the dark energy, there are several 
investigations related to the thawing model for different types of scalar field. All these works are related to individual 
scalar field models and there are not many attempts to study the different thawing scalar field model in a combined 
approach. This is the first such combined approach in python environment. 

This paper is structured as follows. In the section II, we describe the autonomous system of equations for different 
scalar field models and also the method to fix the initial conditions to solve these systems; in section III, we describe 
different cosmological observable that ScalPy calculates for different scalar field models; in section IV, we describe 
different data sets that are used to put observational constraints on different models and the results of the data 
analysis; finally we summarizes our results in the conclusion section. In the Appendix part of the manuscript, we 
briefly describe how to install and to use ScalPy. 


II. DYNAMICS OF COSMOLOGICAL SCALAR FIELDS 

In the present investigation, we confine ourselves to flat FRW universe with scale factor a(t) and given by the line 
element 


ds 2 = — dt 2 + a 2 (t)(dr 2 + r 2 d0 2 ), (1) 

where dQ 2 = (dd 2 + sin 2 (9)dtp 2 ). The scale factor at present is set to be ao = 1. In the late time evolution of 
the universe, the dominant contribution to the energy budget of the universe is from the non-relativistic matter 
(baryons+dark matter) and from the dark energy. In the present paper we are interested in studying various forms of 
scalar field models as the dark energy component in the universe.. 

In following subsection we will look at quintessence, tachyon and Galileon fields and see how we make autonomous 
systems in these frameworks. 


A. Quintessence 

We consider a minimally coupled scalar field and Lagrangian for such a field is given by 

C= l -d^ct>-v(4>), 


( 2 ) 
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where V(cj>) is the given potential. One can calculate the energy momentum tensor and by considering perfect fluid, 
the energy density and pressure of quintessence field is given by 

p .* = + (3) 

P* = (4) 

The equation of motion for scalar field obtained by variation of action with respect to field <j> is 

ij> + 3 Hj> + = 0. (5) 

where the subscript <p denotes the derivative with respect to field <fi. The Hubble equation is given by 


H 2 


1 


where p = p m + p$ is the total density of matter and scalar field. We define the variables x, y, and A as 


i 


( 6 ) 


<t>' 

V6M pl ’ 


A = -M p 


V 


y/V 

V V3HM pl ’ 

r — t /Yii. 
vr 


( 7 ) 

( 8 ) 


Here the prime denote the derivative with respect to the e-folding N (= In a) i.e. <j>' = a(d(j>/da) and subscript ‘0’ 
denote the derivative with respect to field <f>. In terms of these variables 


tip = x 2 +y 2 ,7 = 1 + w = 2 , (9) 

x z + y z 

where is the density parameter for scalar field and w = is the equation of state. We can construct an autonomous 
system using variables x, y and A 


x' = -3x-t- Ay^y 2 + |a;[l + x 2 - y 2 } 

(10) 

y' = ~ x \j\ x v + \v\ l + x2 - y 2 ] 

(11) 

A' = -^A 2 (T- 1)x. 

( 12 ) 

In terms of 7 and A, the above equations becomes '18j, 

i = - 37(2 - 7 ) + A(2- 7 )^ 37 ^, 

(13) 

% = 3(1 - 7 )fy(l-fy), 

(14) 

A' = -n/3A 2 (T- l)^n~ 

(15) 
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B. Tachyons 

Cosmological dynamics for tachyonic field can be obtained by the Dirac-Born-Infeld (DBI) type action fl3j 

Stachyon = - J d^xV1 - (16) 

In natural units, the dimension of the tachyon field is [Mass] -1 . The energy density, pressure and equation of state 
of the tachyon field is given by 


p<t> = 

v(<f>) 

(17) 

\J \ — (j) 2 ' 

p* = 

~V((f>)\Jl — </ 2 , 

(18) 

w = 

-(1 - 0 2 )- 

(19) 


The equation of motion for Tachyon field is 


1 + — ( j ) 2 ) + -^-(1 “ ( f 1 )- 


Again we define following dimensionless variables 

x = H<j/ , y = 


w 


A = —Mr, 


Vs 


VZHM p i 

V„ 


r = v 


(fxf, 


pl 1/3/2 ' x • 

In terms of these variables, equation of state and density parameter f2^ for the tachyon field is 


7 = (i + wj,) = x 2 , 0 . 4 , = 

We construct the autonomous system for tachyon field [a 


y 


\/l — X 2 


l ' = -67(1 - 7) + 2 v /37fi^A(l - 7) 5/4 , 

n^ = 3^(i- 7 )(i-n^), 

V = - V3^VA 2 (l - 7 ) 1/4 (T - 3/2),. 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


(24) 

(25) 

(26) 


C. Galileon 

The large scale modification of gravity which involves an effective scalar field 7r can explain late time acceleration 
of the Universe. This field n is called “Galileon” and its Lagrangian respects the shift symmetry in the Minkowski 
background satisfying n —> n + c and d^cj) —> d^n + b^ where c and b M are constants [13] • The action for Galileon field 
up to third order is jl5| 


S = J d A xV^[^R-\(V „)*(l + JLnn)-V(n)' 


( 27 ) 
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where M p \ = yj is reduced Planck Mass, a is the dimensionless constant and setting it equal to zero givers the 
standard quintessence action. M is the constant with mass dimension and it can be fixed to M p \ by redefining a. 
Again by defining following dimensionless variables, 


x = 


V6HM pl ’ V V3HM p i ’ 
K 
V 


W 


e = 


— 6 — t-Htt . 


A = —M, 


pi- 


r = v 


V-TTK 

w 


we construct the autonomous system 


(28) 

(29) 


3a ; 3 (2 + 5e + e 2 ) — 3a; (2 — e + y 2 (2 + 3e)) + 2\/6y 2 X — \/6x 2 y 2 eX 

4 + 4e + a; 2 e 2 


, y (12 (—1 + y 2 ) (1 + e) — 6a; 2 (2 + 4e + e 2 ) + v / 6a; 3 e 2 A + 2v / 6a; (2 + (2 + y 2 ) e) A) 
V = 8 + 8e + 2a; 2 e 2 

, e (—3a; (—3 + y 2 ) (2 + e) + 3a; 3 (2 + 3e + e 2 ) — 2\f&y 2 X — \/&x 2 y 2 eX ) 

x (4 + 4e + x 2 e 2 ) 

X' = V6xX 2 (1 — T),. 


For e = 0 we recover the autonomous system for the standard quintessence m-m- 
The density parameter f2 m and equation of state for the tt field is given by: 

fim = 1 - z 2 (l + e) - y 2 , 

— 12y 2 (l + e) + 3a; 2 (4 + 8e + e 2 ) — 2\[§xy 2 tX 
= 3 (4 + 4e + x 2 e 2 ) (y 2 + x 2 (l + e)) ' 


(30) 

( 31 ) 

(32) 

(33) 

(34) 

(35) 


D. Solving the dynamical systems 

To solve a cosmological dynamical system, we have to choose a point in time where we set the appropriate initial 
conditions and then let the system evolve till present day. We give initial conditions at around decoupling (a ~ 10 -3 ). 

In case of quintessence and tachyons, we solve the autonomous system of equations in terms of 7 , and A given 
by equations CdD-CESD and (I24l) - (l26|) for quintessence and tachyon respectively. We have to provide initial conditions 
for 7 j, and A i where “i” represents the initial value. As we are considering the thawing class of models, the scalar 
field is initially nearly frozen due to large Hubble friction and hence we choose 7 * to be close to zero. Throughout rest 
of the paper, we fix 7 * = 1(W 4 . A * is the initial slope of the potential and A^ = 0 gives an exact cosmological constant 
behavior. So values for A i different from zero, gives the deviation from cosmological constant. In our subsequent 
calculations, we choose A, to be a parameter of the model. As the contribution of dark energy to be negligible in the 
early time, f2^ should be close to zero initially. But it has to be fixed initially in such a way that for a given 7 * and A 
the system evolves to present day with the appropriate value of where “0” stands for present day value. We check 
that around 10 ~ 9 gives suitable values for So in our calculations, we choose to be the second parameter 
and varies it between 1 x 10 -9 to 10 x 10 -9 . 

For Galileon field, we solve autonomous systems in terms of variables x, y, e and A given by the equations (I30D - (I33D 
and hence we have to provide four initial conditions for Xi, yt, ti and A i respectively. As here also, we consider 
the thawing behavior only, the Galileon field is initially nearly frozen and hence Xi ~ 0. We set a= 10“' in our 
calculations. For to be the same as in quintessence and tachyon, ei and yi are now related by equation (ITU) . So 
we vary €i and as two parameters and tji is fixed subsequently by equation (15H) . We also consider A i as our third 
parameter as it gives the deviation from cosmological constant. So in Galileon models, we have three parameters, 
and A^. 
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FIG. 1: Equation of state w as a function of redshift z for quintessence field (top), tachyon field (bottom left) and Galileon field 
(bottom right). Solid lines corresponds to exponential potential while dashed, dotted-dashed and dotted lines corresponds to 
power law potential with n = 1, 2 and —1 respectively. 


For power law and exponential potential, T is a constant while for other potentials, T evolve as a function of A in 
equations (US,® and ®. 


III. COSMOLOGICAL OBSERVABLES 

Following David Hogg [20] we define observable quantities such as Hubble distance, luminosity distance, angular 
diameter distance and lookback time. The Hubble distance is defined as Dfj = where c is the speed of light and 
Hq is the Hubble constant. Since we have assumed a flat FLRW metric, the line of sight comoving distance is defined 
as 


r dz' 

Hq J 0 h{z') ’ 


D, 


(36) 
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where h(z) = H(z)/Hq is normalized Hubble parameter. The angular diameter distance D A and luminosity distance 
Dl are related to comoving distance as 

D L = (l + z)D m = (l + z) 2 D A , (37) 

where D m is transverse comoving distance which is equal to D c for the flat Universe. We also define lookback time 
to an object as the difference between age of the Universe and the time when the light was emitted from that object: 

tL = ir 0 l (i + z')h(z'y (38) 

The age of the Universe is given by, 

t = 

We also define Om diagnostic [U in our package : 


dz' 


Hq J 0 (1 + z')h(z ’)' 


(39) 


Om(x) = 1 ; x = l + z. (40) 

x 6 — 1 

The time delay distance used in gravitational lensing cosmography is also included in the package. It is defined as j22 ] 

Dm = (1 + Zd) ^ S ■ (41) 

Where Dd is the angular diameter distance to the deflector or lens at redshift Zd and D s is the angular diameter 
distance to the source. Dd s is the angular diameter distance between source and deflector. 

In figure 1 we plotted equation of state, roasa function of redshift for each type of scalar field (quintessence, tachyon 
and Galileon) with power law potential (n = 1, 2, —1) and for exponential potential. The initial conditions we choose 
for quintessence and tachyon fields are (fi^Ai) ~ (2.5 x 10 _ 9 ,0.9) and for Galileon field we choose (fl^e*, Aj) ~ 
(2.5 x 10 -9 , 2.0,0.9). The value of 7 j for all three types of field is taken to be 10 -4 . 


A. Growth of inhomogeneities 

We consider perturbation in matter part only and scalar field is homogeneous. We are interested in structure 
formation to scale much smaller than the Hubble radius so the Newtonian approximation is valid. The density 
contrast of matter field S m = follows the equation at sub-horizon approximation 


5m + 3 H5 m + 4 nGp m 6 m = 0. (42) 

Taking Fourier transform of above equation and defining the linear growth function D and the linear growth rate / 


as 




/ 


D(a)6r, 
d\nD 
din a 


We calculate the linear matter power spectrum defined as 


(43) 

(44) 


P lin (z,k) = A 0 k n °T 2 (k)D 2 n (z). 


( 45 ) 









Redshift z 




Redshift z Redshift z 

FIG. 2: fas redshift z for scalar field (top), tachyon field (bottom left) and Galileon field (bottom right). Solid lines corresponds 
to exponential potential while dashed, dotted-dashed and dotted lines corresponds to power law potential with n = 1, 2 and — 1 
respectively.. 


Here Aq is the normalization constant, n s is spectral index for the primordial density fluctuations generated through 
inflation, D n (z) is growth function normalized such as it is equal to unity at z = 0 i.e. D n (z) = and T(k) is the 
transfer function as prescribed by Eisenstein and Hu [24j] for a mixture of CDM and baryons: 

T(fc)= f^ L ) rb(fc)+ f?K) Tc(fc) , (46) 

\“m 0/ \“m 0/ 

where fl c o is density parameter for CDM and fl m o = fl c o + The form for T),(fc) and T c (k) are given by Eisenstein 
and Hu [241 ]. An important quantity is defined as the fluctuation of mass within the boxes of 8 h~ x Mpc when we 
move from place to place in present day universe.: 



A 2 (k,a) W 2 (k,R)d In k, 


(47) 








9 





^mO 

FIG. 3: la, 2a and 3<r contours obtained for JLA + BAO + Shift(fromCMB) for quintessence field (top), tachyon (bottom 
left) and Galileon field (bottom right) with linear potential. 


where the window function W(k,R) is defined as W(k,R) = 3 ( S ™kR)^ C °(kr ) an< ^ f° r as ' one P u ^ s R = Sh 1 
Mpc. 

Using definition of erg and window function, we fix the normalization constant Aq as [23j 

A 0 = 2tt 2 ctI (^J k ns+2 T 2 (k)W 2 (k,R)dk S j , (48) 

where we take R = 8 ft , -1 Mpc. Fig 2 shows the behavior of fas as a function of z for quintessence, tachyons and 
Galileon fields with exponential as well as power law potential. We used the same initial conditions as used for plotting 
equation of state in figure 1. 


IV. DATA ANALYSIS 


We used our ScalPy package and MCMC hammer ‘emcee’ da to constrained model parameters of quintessence, 
tachyon and Galileon fields as well as wCDM, wo-u^CDM and GCG models(See Appendix). We used following data 
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FIG. 4: 1 < 7 , 2a and 3a contours obtained for JLA + BAO + CMB Shift (dashed) and JLA + BAO + CMB Shift + fas (solid) 
on model parameters for wow a CDM model. 


sets for our analysis: 

• 31 binned distance modules fitted to JLA sample given by by M. Betoule et al. J2f| 

• Combined BAO/CMB constraints on the angular scales of the BAO oscillations in the matter power spectrum 
measured by SDSS survey, 6dF Galaxy survey and the Wiggle-z survey. We used covariance matrix for this 
provided by Giostri et al [26| . 

• We also used measurement of shift parameter of CMB as given in Shafer et al [27| ■ 

• Measurements for fas by various Galaxy surveys like 2dF, SDSS, 6dF, BOSS and Wiggle-Z. A compilation of 
different measurements for fa§ has been provided by Basilakos et al. |28| 

The results for data analysis for quintessence, tachyon field and Galileon field with linear potential are shown in figure 
(3). We used h = 0.7 in our analysis. We show the la, 2a and 3a contours in — w plane. The ACDM behavior 
(w = —1) is allowed for all these types of scalar fields within la region with considerable range of fi m o values. One 
can also see that the constraints on present day equation of state wq and density parameter for matter £l m o are same 
for all three type of scalar fields with linear potential. This result also holds for other power law potentials. 

In figure (4) and (5) we show the constraints on model parameters for woU'aCDM and wCDM models respectively. 
The concordance ACDM model is within la error bar when we combine all the four data sets used. The constraints 
on the model parameters of GCG model is shown in figure (6). The parameter A s is related to present day equation 
of state wq as wq = — A s . We see that wq = — 1 is again allowed within la region for GCG as well. 


V. CONCLUSION 

To summarize, we present a python package ScalPy to study the late time cosmology with different kind of scalar 
field models. To begin with, we consider only the minimally coupled scalar field and a flat FRW universe. We also 
neglect the radiation component as it will not affect the results in the late time. But one can trivially incorporate 
radiation component in ScalPy. As thawing models of scalar fields are more consistent with different observational 
data, we confine ourselves only to the thawing class of fields where the scalar held is initially frozen due to large 
Hubble friction and starts very close to w = —1. ScalPy solves the autonomous system of equations for different 
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FIG. 5: la, 2a and 3<r contours obtained for JLA + BAO(dotted),JLA + BAO + CMB Shift (dashed) and JLA + BAO + 
CMB Shift + fas (solid) on model parameters for wCDM model. 




FIG. 6: la, 2a and 3cr contours obtained for JLA + BAO + CMB Shift (solid) and JLA + BAO + CM B Shift + fas (dashed) 
on model parameters for GCG model. Present day equation of state wo is equal to — A s in GCG model. 


kind of scalar field models with power-law and exponential potentials. It calculates various observables related to 
background expansion as well as to the growth of structures in universe. In this paper, we consider only the linear 
regime in the sub-horizon scales where effect of dark energy perturbation is negligible. Finally we integrate it with 
MCMC sampler “emcee” to constrain various scalar field models with presently available data. 

To conclude, we present a python package to study various scalar field dark energy models. This is probably the 
first package exclusively for the scalar field case. At present, we confine only to the canonical scalar fields and flat 
FRW models. But in near future, we shall incorporate various other scalar field models like coupled quintessence, 
scalar tensor models as well as non-flat FRW models. We shall also include the observables related to super-horizon 
scales where dark energy perturbation can play an important role. 
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Appendix A: Installation 

This package is written in python and it solves dynamical system for three types of minimally coupled scalar 
fields- quintessence, tachyon and Galileon with exponential and power law potential. Given the initial conditions, this 
module solves for cosmological observables such as Hubble parameter, luminosity distance, angular diameter distance, 
lookback time, growth function, growth rate, fag(z) and power spectrum (with Eisenstein-Hu transfer function as 
well as BBKS transfer function). We have also defined all above cosmological observables for standard cosmological 
fluids such as ACDM, wCDM, wow Q CDM and generalized Chaplygin gas (GCG) H^]. Table 1 describes different 
cosmological models we have used in the code. 


TABLE I: Different parametrizations for dark energy models 


model 

parameters 

Hubble equation 

ACDM 

wCDM 

wowaCDM 

GCG 

^mO 
^mO) IV 

wo , w a and Q m o 
fimo, A s and a 

= ' 2 ') 3 “1“ (1 — ^mo) 

= Vn m0 (i + *) 3 + (i - n m0 )(i + z)- 3 ^" 0 

= \/n m o(l + z) 3 + (1 - Hmo)(l + z) 3( - 1+w 0 +™<d exp(—3w a z/(l + z)) 
= V«™o(l + z) 3 + (1 - n m0 )(A s + (1 - A s )(l + Z )3(1 +o))1 /(1+ Q ) 


To install the code, one can simply use pip or download the full source code and run ‘setup.py’ file. On a linux 
system where pip is installed, scalpy can be installed as 
user@computer > sudo pip install scalpy 

To install it from source code which is downloadable as .zip file from github repository, simply download the source 
from 

https://github.com/sum33it/scalpy 

After extracting the file, move to the folder ‘scalpy-master’ and run ‘setup.py’ file 
user@computer > unzip scalpy-master.zip 
user@computer > cd scalpy-master 
user@computer > python setup.py install 


Appendix B: Structure 

The package consists of following main modules linked with each other: 

1) scalar.py This module solves the dynamical system for scalar fields such as minimally coupled quintessence, 
tachyon and Galileon. We have included two types of potential for each scalar field: power law U(^>) = Kef/ 1 and 
exponential V(cf>) = exp (a<j)) where n,a and K are real constants 

2) fluids.py This gives different observables for standard cosmological models such as ACDM , wCDM , wow a CDM 
and GCG In this hie also, different observables such as luminosity distance, Hubble parameter, angular diameter 
distance, growth rate, growth function, power spectrum are defined 

3) solver.py For a given values of wq and o, this module solves for the initial conditions needed at decoupling for 
scalar held models. 

4) transfer_func.py In this module, we have defined transfer functions given by Eisenstein and Hu as well as by BBKS. 
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Detailed online documentation with examples can be found at following url: 
http://ctp-jamia.res.in/download/scalpy 

One can also read readme.rst file for quick Introduction provided in tarball downloaded from github repository. 
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